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p 

Let  S  :=  if  be  the  space  of  bivariate  piecewise  polynomial  functions 

JC  /  u 

in  <?,  of  degree  <  k,  on  the  mesh  A  obtained  from  a  uniform  square  mesh 
by  drawing  in  the  same  diagonal  in  each  square. 

de  Boor  and  HOllig  have  given  the  following  upper  bound 

m  <  m(k)  :=  min{2(k-p),  k+1} 
for  the  approximation  order  m  of  S. 

In  this  paper,  the  lower  bound 

m  >  m(k)  -  2 

is  demonstrated.  This  result  is  close  to  de  Boor  and  HOllig's  conjecture 
that  m  never  differs  from  m(k)  by  more  than  1. 

Incidentally,  the  approximation  order  of  irj  .  is  shown  to  be  4. 
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SIGNIFICANCE  AND  EXPLANATION 


Univariate  splines  have  been  proved  quite  useful  in  practice.  However, 
if  one  wants  to  fit  a  surface,  or  solve  a  partial  differential  equation 
msnerically,  one  would  naturally  think  of  using  multivariate  splines.  Here 
splines  still  mean  piecewise  polynomial  functions.  In  this  respect,  a  basic 

question  is  to  ascertain,  for  a  given  mesh  and  a  family  S  of  splines  on 

*1 

A',  what  its  optimal  approximation  order  is.  This  question  is  challenging 

-/UX.  fx-J) 

even  for  a  regular  triangular  mesh  ^^as  soon  as  one  demands  that  the 

/ 

approximating  functions  have  a  certain  amount  of  smoothness.  The  report 
records  a  step  toward  answering  the  above  question.  - 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
8Uftnary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


APPROXIMATION  BY  SMOOTH  BIVARIATE  SPLINES 
ON  A  THREE-DIRECTION  MESH 

Rong-qing  Jia 

1 .  Introduction 

In  this  paper  we  study  approximation  order  of  smooth  bivariate  splines  on 
a  three-direction  mesh.  The  work  in  this  respect  was  initiated  by  de  Boor  and 
DeVore  IBD]  and  de  Boor  and  Hdllig  [BH  1,2,3].  Here  we  follow  them  and 
introduce  some  notations.  Let 

A  U  {x  e  R2;  x(1)  =  n,  x(2)  =  n,  or  x(2)-x(1)  =  n>  . 
nes 

Namely,  the  mesh  A  is  obtained  from  a  uniform  square  mesh  by  drawing  in  the 
same  diagonal  in  each  square.  Let 

o  P 

S  :=*  w,  s»  ir.  .  n  C 
k, A  k, A 

be  the  space  of  bivariate  pp  (piecewise  polynomial)  functions  in  CP,  of 
total  degree  <  k,  on  the  mesh  A.  Also,  by  we  denote  the  space  of 

polynomials  of  total  degree  <  k.  we  are  interested  in  the  approximation 
order  of  S.  The  approximation  order  of  S  is,  by  definition,  the  integer 
m  for  which  the  following  hold',:  For  all  sufficiently  smooth  function  f. 


dist(f ,Sh)  -  0(hm) 


while,  for  some  C  -function  f, 


dist (f , S^)  ?  o(hm)  . 

Here,  the  scale  (S^)  of  approximating  spaces  is  generated  from  S  by  simple 
scaling. 


s.  j*  a.  (S) 
n  n 
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(0.  f)(x)  :=  f(x/h),  all  f,  x,  h 
n 

de  Boor  and  DeVore  have  given  the  following  lower  bound  for  m  (see 
IBD] ) s 

m  >  p+2  in  case  p  <  p(k)  s=  [_(2k-2)/3_J 
In  contrast,  S  has  approximation  order  0  for  p  >  p(k). 

An  upper  bound  for  m  has  been  obtained  by  de  Boor  and  HOllig  (see  [BH 
3;  Theorem  3] : 

m  <  m(k)  :=  min{2(k-p),  k+1} 

de  Boor  and  HOllig  also  show  that  the  approximation  order  of  n*  is  3 

& 

rather  than  4  (see  [BH  2]).  Thus  the  approximation  order  of  S  may  differ 
from  m(k)  by  1.  Based  on  those  investigations,  de  Boor  and  HOllig  raised 
the  following 

Conjecture  ( [BH  3] ) .  The  approximation  order  of  S  =  w?  never  differs 

A 

from  its  upper  bound  m(k)  by  more  than  1. 

In  this  paper,  we  shall  show  that  the  approximation  order  of  S  = 

A  f  A 

never  differs  from  m(k)  by  more  than  2.  The  proof  of  this  result  will  be 

based  on  a  quasi-interpolant  scheme.  For  the  record  we  state  the  following 

Theorem  1 .  Suppose  that  B  e  S  with  supp  B  finite.  If  the  map 

T  :  p  t*  l  p(  j)B(  -j) 

ze*2 

is  one-to-one  and  onto  *  ,  then 
-  n  - 

dist(f,Sh)  =  0(hn+1 ) 

for  all  sufficiently  smooth  functions  f . 

The  argument  in  [BH  1;  Section  6]  essentially  gives  the  proof  for  Theorem 
1.  We  do  not  need  to  repeat  the  proof  here. 

To  construct  an  element  B  B  S  with  the  property  required  by  Theorem  1, 
we  shall  emply  box  splines,  which  were  introduced  by  [BD]  and  [BH  1].  In 


section  2,  we  develop  some  preliminary  results  from  univariate  B-spline 


theory.  In  section  3,  we  elaborate  some  properties  of  box  splines  on  the 

three-direction  mesh  A.  In  section  4,  we  construct  an  element  B  e  S  with 

the  property  required  by  Theorem  1,  and  therefore  prove  our  main  results.  I 

section  5,  we  show  that  the  approximation  order  of  ^  is  4.  This 

0 

illustrates  that  the  approximation  order  of  *  might  be  exactly  k  when 

K/  A 

k  -  2p+2. 


2.  Some  preliminary  results  from  univariate  B-spllne  theory. 

m 

Let  t  -  (t . )  be  a  knot  sequence.  Recall  that 

*  1  "* 

Mi  k(x>  *"  fclt4#...#t1+k]  (— x)J  1 

is  a  normalized  B-spline  of  order  k  for  each  i  e  Z.  Also  we  write 

"i,*00  ui+k  *  ti)"i.k,,t)^  • 

If  p  is  a  polynomial  of  degree  <  k,  then 

P  *  l  , 

iez  1  1 #x 

where  is  the  linear  functional  defined  by 

A.  kf  -  l  C-,k-,-j^1-j,f(j,CT1, 
i,k  j<k 


♦i,k(x)  s“  (ti+rx)  *•*  (ti+k-rx)/(k*1,‘ 

and  e  (t^,t^+k)  (see  IBF] j  also  [B]).  Now  suppose  *  i,  all  i  e 
Then  “  Mj,  fc  end 

♦  .  .  (x)  -  (i+1-x)  • • • ( i +k- 1 -x ) / ( k- 1 ) I  . 

1  #IC 

l  It  is  easily  seen  that  there  exist  unique  constants  a^  ( 1*0, 1, . . . #k-2) 
such  that 

k-2 

♦A  Jx>  “  l  *•*  (*+k-2-x)/(k-2) ! 


Comparing  tha  coefficient  of  x 


on  both  sides  of  (1),  we  obtain 


If  f  is  a  polynomial  of  degree  <  k-2,  then 


P  *"  I  fd)N  .. 

iex  i,k  1 


is  also  a  polynomial  of  degree  <  k-2.  On  the  one  hand. 


f(i)  -  X 


i,k-1p  ' 


On  the  other  hand. 


P  *  l  kP)N.  .  • 

iex  i,K 

Pick  T  6  (i+k-2, i+k-1 )  and  calculate  X  p  as  follows: 
i  1  ,  K 

1,K  l,k  i  1 


l  (-) )(k‘2“j)(r 
k-1  i,k 

-  I  (-)k’1“3t-Yaik.1^>(xi,p^>(Ti 

j<k-1  4-0  *'k  *  1 

ky2  a  y  (->k-2-Vk’2_3)<T  )n(j)(T  > 

,io  *'k-’  J-1  ‘  ’  <v 


)] 


k-2 

I  -tk. 

4-0 


1  Xi+4,k-1p 


Therefore 

L  f(1)N1,k-i  "  p  =  (Xi,kp)Ni,k  "  a4,k-1  ,I.(Xi+4,k-1 


iex 


iex 


iex 


X  L  f(1+t)Ni,k  -  A  a*.k-i  J  f(1)Ni,k(*+A) 


iex 


4-0 


4-0 

We  have  proved 

Lemma  1 .  For  any  polynomial  of  degree  <  k-2 

k-2 


iex 


l  k_,  -  l  *<i)(  l  V-1N,  k<-+*>> 

iex  i,k  1  iex  4-0  *'k  1  i,k 


p)i 


3.  Box  aolines  on  a  three-direction  mesh. 


As  defined  in  [BH  1],  the  box  spline  M„  is  the  distribution  on  rf" 


given  by  the  rule: 


M *{  l  X(i)5.)dX 
S  [0,1]*  i-1 


for  some  sequence  3  :*  in  •  1°  our  case*  m  *  2.  I*t  be  the 

unit  vector  along  x^-axis  (i*1/2),  and 

:■  ey  d2  :*  e^+ej*  d3  :*  e2  • 


.  r+s+t 


For  positive  integers  r,  a  and  t*  let  3  *  (C.)^  be  the  sequence  in 

K2  given  by 


«r  *  V  «r*1 


5r+s  "  d2  and  ^r+s+1 


^r+s+t  *  d3 


From  now  we  will  write  ML  .  «.  instead  of  M„.  Caution  1  Our  notation  is 

r/B# w  s 

slightly  different  from  [BH  3J .  In  [BH  3],  d2  *  e2  and  d3  *  ei+e2*  Thus 
our  H,  _  „  is  just  M  .  in  the  sense  given  by  [BH  3} . 

Hie  smoothness  of  M_  _  ..  depends  on  the  direction  multiplicities. 

*#■#  V 


From  [BH  3]  we  have 


M  „  e  i£«  c  c<d-1) 
r,S,t  •  ~ 


with  d  •  min{r+s,s+t,t+r}  -  1. 


Now  we  define 


r-2  0  8-2  0 


Br  «  t<xi'x2)  S“  I  l  l  **•  l 

T,a,t  1  1  X  -0  X  -0  u  *0  u 


x,-°  xr-r°  v°  vr°  vr°  vr° 


(a.  •••a.  a  •••a  .a  ...a  x 

V  Vi'**’  V1  Vi'**'  V’  Vi'1*’ 


for  (r,s,t)  with  min{r,s,t)  >  1  , 


where  a  has  the  meaning  determined  by  ( 1 ) . 


-5- 


reason  for  introducing  B_  „  *.  will  be  clear  after  we  prove 

r  tQfZ. 


following 


Lemma  2.  For  any  bivariate  polynomial  of  degree  <  r+s+t-2,  we  have 

1°  I  2  p(j)(Br;S>t  -  »r,.,t_1)(-JH  -  o  . 


2°  Dfo!  I  2  POHBri5it  -  Brj,.1it)(-J)l  *  0  , 


3°  l»,«!l,>;[!]Pli)(Bril,t-.r.lilitl(-)l)-.  • 


Here,  as  usual,  j  **  (j^,J2)  e  z2»  x  *  (x1,x2>  e  R2,  D±  - 
VAf  -  f-fC-e^  (i  -  1,2). 

Proof.  By  symmetry,  2°  and  3°  follow  from  1°.  Thus  we  only  need  to 


prove  1  .  Suppose 


J  ^  .  *  I  b.M  _  .(•+!)  . 

r ,s,t-1  L  _  i  r,s,t-1 

ie*2 


i,  by  the  definition  of  B_  „  .  ,  we  have 

r ,  a ,  c 


*  ^  “  I  b.  I  a.  ,M  (B+i+Jle..)  . 

r,8,t  u  i  ,  l,t-1  r,s,t  2 


.2  1-0 


Hence 


ie* 


[li  !,  l^W’pHlXt  I  *l,t-,’Vo.t<“1*lt*V1  -  "o.o.t-i'-i*1” 

1  .  W  w 


For  any  test  function  $,  one  can  easily  check  that 

<"0.0.t’*>  "  ^  Mt(x2,+(0'x2)dx2  * 


t-2 


<  I  I'llV’/pIDIl!  •,,t-i%.0,t,-jtit“21  -  "o.o.t-l1"1*1”'  *y 

j etr 


t-2 


/  I  (V1J(71+V2)8p)(j)[  I  ai/t.1Mt(x2-j2+i2+£)-Mt_1(x2-j2+i2)]*(0,x2)dx2  =  0 

je*2  1-0  ' 


by  Lemma  1,  since  V^(7^+72>8p  has  degree  <  t-2.  This  completes  the  proof 
of  Lesssa  2. 


4.  Quasi-lnterpolant  scheme. 

In  the  following,  r,  s  and  t  are  always  integers.  Let 
I  s»  {(r,s,t)  |  r+s+t  -  2p+4  and  2  <  r,s,t  <  p+1} 

{(r,s,t)  |  r+s+t  »  2p+3  and  2  <  r,s,t  <  p+1} 

J2  s-  {(r,s,t)  |  r+s+t  *  2p+3  and  2  <  r,s,t  <  p} 

K  {(r,s,t)  |  r+s+t  -  2p+2  and  2  <  r,s,t  <  p}  . 

We  have 

I  -  { (r,2,t) |r+2+t-2p+4,2<r,t«p+l}  U  {(r,s,t) |r+s+t»2p+4,3<s<p+1,2<r,t<p+1 } 

-  {(P+1, 2, p+1)}  u  {(r,s,t) |r+s+t-2p+4,3<s<p+1,2<r,t<p}  u 
{ (r,s,t) |r+s+t-2p+4,3<s<p+1,2<r,t  and  max{r,t}»p+1 } 

-  {(p+1, 2, p+1)}  u  {(r,s,t) |r+s+t-2p+4,4<s<p+1,2<r,t<p}  u 

{ (r,s,t) |r+s+t-2p+4,3<s<p+1,2<r,t  and  max{r,t}-p+l}  .  (4) 

Similarly, 

J1  -  {(r,s,t) |r+s+t-2p+3,  3<s<p+1  and  2<r,t<p}  u 

(5) 

{ (r,s,t) |r+s+t-2p+3,  2<s<p,  2<r,t  and  max{r,t}  -  p+1}  , 

and 

J2  «  { (r,s,t) |r+s+t  -  2p+3,  3  <  s  <  p  and  2  <  r,t  <  p}  .  (6) 

Therefore 

U|  ;j2|  -  Ij,!  +  |k|  -  i  .  (7) 


Here,  by  |B|  we  mean  the  cardinality  of  the  set  E 


In  the  following  we  use  the  convention  that  the  empty  sum  has  value  0. 
Now  we  construct  B  as  follows: 


BS=X  £  B  *  1  B 

/  *#s»t  .  r,s,t 

(r,s,t)ei  (r,s,t)eJ1 


l  B  +  l  B 

r»S,t  .  t  T,B,t 


(r,s,t)eJ 


(r,s,t)eK 


(8) 


i 


Lemma  3.  £  B(  -j)  =  1. 

je*2 


Proof.  From  [BH  1]  we  have 


Then  (2)  and  (3)  yield 


therefore 


E  'j)  *  1  ' 

je* 


l  'I'  *  1  • 

je* 


l  B(  -j)  *  III  -  |J  I  -  I J- |  +  |K|  -  1  . 

je*2 

The  following  lemma  plays  an  essential  role  in  this  paper. 

O 

Lemma  4.  For  k  *  2p+2,  B  e  «[  .  and 

~  K  f  A  ~ 

p  -  £  p( j )B(  -j)  is  a  polynomial  of  degree  <  deg  p 

je*2 

for  any  polynomial  of  degree  <  k- 1 . 

Proof.  We  first  show  that 


^1  *^2  r  2 

Di  Dj  [£  p(  j ) B(  -j)]  is  a  constant  for  any  e  * 


je* 


with  >  0,  g^  >  0  and  =  de9  p  *  k-1 


(9) 


-8- 


Fix  and  q2»  Consider  the  following  index  sets: 


E1  :*  {(r,s,t)  |  r  >  q1  and 

E2  :  =  {(r,s,t)  |  r  <  <?1  an<i 

E3  :*  {(r,s,t)  |  r  <  q1  and 

E.  :*  {(r,s,t)  |  r  >  q.  and 

4  1 

Then  {ea  j  i  *  1, 2/3,4}  forms  a  partition  of 
sufficient  to  show  that 


t  >  q2> 
t  <  q2> 
t  >  q2> 
t  <  q2>  • 

Z^.  To  prove  (10)  it  is 


^2  v  r  r  r 

D,  D  l  p<j>(  l  B  -  l  Br 

je*2  (r,s,t)einEi  '  '  (r,s,t)eJ1nEi  ' 


2.  B  + 

(r.B.tleJ  "I  r'''  (t,s,t)exnE 


I  Br.t)(,'j) 


is  a  constant  for  each  i  -  1,2,3  or  4.  Thus  we  have  to  split  our 
consideration  into  the  four  cases:  i  *  1,2,3  or  4. 

Case  1*1.  Then  r  >  q.,,  and  t  >  q2.  We  have 


q1  q2 


qi„q2 


D1  °2  \  P<*)Mr,s,t(-:»  *  \  {71  V2  P,(^Mr-q1,s,W2(^)  ' 

ql  q2 

which  is  a  constant,  since  V,  p  is  a  constant.  It  follows  that 
q1  q2  r 

D.  D_  )  p(j)B  .  (■•-j)  is  a  constant  for  any  (r,s,t) 

1  2  *  r,s,t 

with  r  >  q^  and  t  >  q2 


Hence 


D,V  I  P(j>  <  l  Br  s  t  "  ^  Br  s  t 

jez2  (r,s,t)einE1  '  '  (r,s,t)eJ1nE1 


i  Br  s  t  +  i  Br  s  t)('“j) 

(r,s,t)eJ2nE1  ,B,Z  (r,s,t)eKnE1  T,a,z 


is  a  constant 


Note  that 


Case  i  *  2.  In  this  case,  r  <  a,  and  t  <  q„ 

1  d 

(p+1,2,p+1)  f.  E2  . 


This  is  true,  because  r  <  q1  and  t  <  imply  that 


r+t  <  q1+q2  <  2p+1 


Now  (4),  (5)  and  (6)  tell  us  that 


r,s,t 


r,s,t 


(r,s,t)einE„  (r,s,t)eJ  nE  (r,s,t)eJ,nE 


B  .  + 
r,s,t 


12 


2  2 


P+1 


(r,s,t)eKnE 


l  Br  s  t  "  I  l  (Er  s  t  “  Br  s-1  t 

i  r/S,c  _ .  r,s,t  r,s  i,t 


)  + 


S=4  2<r,t<p 

r+t=2p+4-s 
r<q1,t<q2 


P+1 

I  l 

s=3  2<r 

max{r, t}=p+1 
r+t=2p+4-s 


P+1 


,t  ( Br , s , t ”Br , s - 1 , t }  ‘  J3  2<rJt<p  (Br,s,t*Br,s-1,t 


r+t=2  p+3-s 


r<q.,/t«q2 


r<q1#t<q2 


while 


qi  q2 


D,  °2  !  I  P(JXBr>I>t-»r<^1(t)(-J))  - 


jez 


D?  'DI2  VD2[  l  2 


jez 


by  Lemma  2.  Therefore 


qi_q2. 


D,  °2  f  t  I  Br  s  t  '  ^  Br  s  t 

jez2  (r,s,t)exnz2  r,8,t  (r,s,t)eJ1nE2  r,s,t 


l  B  + 

(r,s,t)eJ2nE2  r'8'  (r,s,t)eKnE2 


l  B  s  t)(*_j)]  = 


Case  1=3.  Then  r  4  q_  and  t  >  q 


We  have 


VD2  -D1D1  °2  -VVVV  °2 


qrr  q,-r-t  I  q  -r-t  q  -r  q 

1  (  *  K 


*  'O.  4  *  ^  « 
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D*[  l  (-D  (  4  )(D  +D  )*b  1  2  ]  , 

4=q1+q2-r-t+1 


(11) 


where  the  differential  operators  H  „  and  G_  .  are  defined  as  follows 

IT  #  S  Fit 

q  -r  m 

1  q.-r-4  q  -r  , 

"t,.  •’  J  (-”  (  1  )‘W 


4=8 

q^q^r-t 


l-s  Vq2_r_1 
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r,t 


q  -r-4  q-r 

l  <“1>  (  *  2) 

4=0 


D2 


For  the  third  term  in  (11)  we  observe  that 


r  4  Vq2"r“*  n 

VW^  1  I 

jear 


,  q,*q2-r-‘ 


I  •’’[‘VV  V  p)(j)H 


je* 


0,s-A,t-(q1+q2 


-r-4) ( 


is  a  constant  for  4  e  tq1+q2”r"t+l*B”'U  •  because 


s-4  >  0  and  t  -  (q  +q  -r-4)  >  0 


Thus  we  can  omit  the  third  term  of  (11)  in  the  following  discussion.  Now  we 


want  to  prove  that 


l  P<j><  l  -  l  -  l  ♦  i  IIH  D'(D ,«  )S  *  Or,t<OjJBr,e>tC-i) 

1  Tne*  T  n  P  "T  Ap  Va  W  9  w  w  w 
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is  a  constant.  Let 


®  *-  l  P<j><  I"  I  -  I  +1  )(Hr  sD1(D1+D2)3Br  s  t(*'j)) 

O  T  T  T71  r»S  1  1  *  I,8,Z 
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inE3  J1nE3  J2nE3  KnE3 


V  :=  l  p(j)(  I  -  l  -  l  +  1  )(Gr  tD1D2)Br  s  t(“j)) 


jez2  InE3  JinE3  J2nE3  KnE3 


Then  we  can  argue  separately  for  U  and  V.  Interchange  s  and  t  in  (4) 
(5)  and  (6).  We  can  write  down  (cf.  (10)) 

(-j) 


°  - 1 2  p(j,v..p+.  Dr,(D.*D2>p*’Bp+i,p+i,2' 


jez 


p+i 
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r+s=2p+4-t 

r<q1#t>q2 


P+1 
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r+s=2p+4-t 


r<qipt>q2 


p+i 


+  ^  p(j)  ^  I  Hr  sD1<D1+D2)S(Br  s  t  '  Br  s  t-1)(*'j)  * 

t-3  2<r,s<p  r'5  r'S-‘  t'»'t 


r+s*2p+3-t 

r<q1*t>q2 


However,  Hp+1>p+1  =  0p  because  q1  <  2p+1  <  (p+1)  +  (p+1).  For  the  second 


term  of  the  above  expression,  we  have 


1  2  P(j)  A  D  Hr,aD1(D1+D2>  (Br,s,t  '  Br,s,t-1 > < ^ 

.2  t“4  3sr,ssp 

3-  r+s»2p+4-t 

r<q1,t>q2 

‘  J.  I  l  2P<1,,Br,e.t  ’  ‘  0 

t*4  34r,ssp  ,2 

r+s«2p+4-t  3 

r<q1 ,t>q2 

by  Lemma  2.  The  third  and  fourth  terms  of  (12)  are  also  zero  by  the  same 
argument,  lb  us  U  *  0.  Similarly  we  can  show  V  *  0.  Therefore 


q1  q2 


D,  V  t  1  P<jH<  I 


H1E3  J1nE3  J2nE3  K0E3 


l  *  l 


is  a  constant. 


Case  i  »  4.  In  this  case  r  >  q^  and  t  4  q2»  and  the  argument  is  as 


in  Case  3. 


So  far  we  have  proved  statement  (9).  Now 

P  -  I  p(j)B(—  j) 


is  a  polynomial  of  degree  <  deg  p.  For  (q3,q2)  with  q^  >  0,  q2  >  0  and 


q1+q2  "  <3eg  p 


we  have 


qi  q2 

■*  'n  4  , 


*1*2 


*1  *2 
n  1  r» 


V1  V(P  ■  2  PCj)B(-j))  -  V/Vj  p  -  l  p(j)(v/v/B(.-j)) 

j 


vi1y22  p  ■  ^  (vr722p)(j)B(,_j)  • 


q1  q2 
n  'n 


*1  *2 


However,  ?2  p  is  a  constant.  Hence 


q1  q2 
rt  1  rt  4 


*1  *2 


l  <VYP)(j,B(-i)  “  7i  Yp 


by  Lemma  3.  Therefore 


V1  V2  (p  -  l  p(j)B(*-j))  -  0,  for  any  (g^q^  with  q1>0,q2>0  and  q^q 

je*2 

This  shows  that  p  -  J  p( j )B { •— j )  is  a  polynomial  of  degree  <  deg  p. 

je*2 

Lemma  4  is  proved. 


Now  we  can  prove 


Theorem  2. 


EJdiEtl 


T  defined 


T  :  p-+  l  p(j)B(—  j)  ,  p  e  s 

jei2 

is  one-to-one  and  onto  ir. 

- — — . . —  k- 1 

Proof,  w.  „  is  a  linear  space  of  finite  dimension,  and  T  is  a  linear 
k-1 

mapping  from  ^  to  (  by  Lemma  4.  If  p  ^  0,  then  deg  p  >  0.  Lemma 

4  tells  us  that  £  p(j)B(*-j)  has  the  same  degree  as  p;  that  is 
j 

X  p( j ) B < ■ — j )  ?  0.  This  shows  that  T  is  one-to-one.  Since  w  is 

2  * 

je* 

finite-dimensional,  T  is  also  onto.  The  proof  of  Theorem  2  is  complete. 

Now  combining  Theorem  1  and  Theorem  2  gives 

P 

Theorem  3.  If  k  -  2p+2  and  S  *  then 

distff^)  -  0(hk) 

for  any  sufficiently  smooth  function  f . 

Remark.  From  the  above  arguments  we  see  that  Theorem  3  remains  true  for 


k  >  2p+2. 

We  show  in  Section  5  that  the  approximation  order  of  ^  is  4.  Thus, 
in  general.  Theorem  3  cannot  be  improved. 

For  the  general  case,  we  also  have 

Theorem  4.  If  S  *  wf  .  and  p  4  p{k)  :»  L(2k-2)/3j,  then 

dlst(f , S^)  -  0(hra(k)“2) 
for  any  sufficiently  smooth  function  f. 


Proof.  Proa  [BH  1]  we  already  know 


0+2 

dlst(f,Sh)  -  OUT  )  . 

If  2k  <  3p+4,  then 

m(k)  -  2  <  2{k-p)  -  2  -  2k  -  2p-2  <  p+2  . 

Hence  Theorem  4  holds  for  2k  <  3p+4.  If  k  >  2 p+2,  then 

a(k)  -  2  <  k-1  . 

thus  Theorem  4  follows  froa  Theorem  3.  How  assume  2k  >  3p+5  and  k  <  2 p+2 
Let 

0  t-  2p+2-k,  k'  tm  k-3o,  p'  :■  p-2o  . 


P  *-“►  l  p( j)B(»-j) 

jez2 

is  one-to-one  and  onto  W2k-2p-3*  Now  ^eorem  ^  gives  the  required  result 

For  any  sufficiently  smooth  function  f, 

dist(f , S  )  -  0(h2k"2p"2)  . 

n 

This  ends  the  proof  of  Theorem  4. 


roximatlon  order  from  bivariate  C'-ouartics 


In  this  section  we  will  show  that  for  S  ■  if.  .  and 

4,  A 

f  s  ►-♦  x2x^  /  x  -  (x1#x2)  e  R2  , 
there  exists  a  positive  constant  such  that 

dist(f,S.  )  >  const»h*  . 

n 

To  this  end  we  shall  follow  [BH  2]  and  discuss  B-nets  in  the  following. 
Given  a  triangle  T  with  vertices  U,  V  and  W,  we  associate  each 
point  x  with  its  barycentric  coordinates,  i.e.  with  (u,v,w)  for  which 

x  ■  uU  +  W  +wW,  and  u+v+w  -  1  . 

Any  polynomial  p  of  degree  <  n  can  be  represented  by 

p '  bi^k 

with 

nl  ilk 

♦ijk'*’  ■■  rnrirr u  v  *  • 

where  bj^  are  uniquely  determined  by  p.  This  representation  gives  rise  to 
a  function 

b  1  xi  jjc  ^ijk'  Xijk  !"  jV+kW)/n  and  i+j+k  *  n 
This  function  is  called  the  B(ernstein  or  ezier)-net  for  p  (with  respect  to 
T).  (See  (BH  2] .) 

To  a  given  function  f  e  ^  we  associate  a  function  b^  so  that  b^ 
is  defined  on 


»-  <R/4)‘ 


6- 


and  bf  agrees  with  the  B-net  for  f  on  each  triangle  cf  A.  Obviously, 

bf  is  well  defined.  We  also  call  bf  the  B-net  for  f  with  respect  to  A. 

0 

Let  us  now  introduce  some  linear  functionals  on  it.  ..  Define 

4,  A 

^(m,n)f  im  b^(m+  i~,n)  +  bf(m+  ^-,n)  -  bf(m+  ~-,n-  ~)  -  bf(m+  ^-,n+  ~)  , 

^^(m,n)f  bf(m,n+  ^~)  +  bf(m,n+  ~)  -  bf(m-  ~,n+  -  bf(m+  ~,n+  ~)  , 

Xi3m,n>f  s“  bf(m+  'V''n+  i4^)  +  bf(m+  4'n+  4°  *  bf(m+  ^'n+  4°  " 

b.(m+  "7,n+  ■“) ,  i  “  1,2, 3,4;  m,n  B  Z 
f  4  4 


Let 


and 


Aij  !=  {Xlj'n)  I  m*n  e  z>  *  i  -  1*2, 3, 4,  j  «  1,2,3 


4  3 

A  U  U  A  . 
i-1  j-1  13 


If  f  e  it 


4, A 


then  Abf  *  0  for  any  A  6  A  (see  [F]  and  [BH  2]) 


We  extend  each  A  B  A  to  the  continuous  linear  functional  Al  on  C(R2) 
with  the  aid  of  the  local  linear  map  I  which  associates  f  with  the  unique 
element  If  of  it^  ^  which  agrees  with  f  on  J4«  Let  T  be  the  mapping 
f  t-»  bJf  for  f  e  C(S?),  and  let  Tj  be  the  shift  operator  f  *— *■  f(»+j). 

We  have  the  following 

Lemma  5.  T  is  a  linear  mapping  and  commutes  with  any  shift  Tj ,  j  B  Z2. 

Proof.  It  is  obvious  that  T  is  a  linear  mapping.  To  prove  the  second 


statement  we  first  show  that  I  commutes  with  any  Tj .  Indeed, 

Tj { If) ( i)  -  If(i+j)  ■  f(i+j)  for  any  i  e  J4  , 

I (Tjf ) (i )  =  I(f(«+j))(i)  -  f ( j+i )  for  any  i  B  J4  . 
This  shows  that  Tjl  *  ITj.  N®xt,  have  to  show  that  the  mapping 

bg  ,  9  e  .“  4 

commutes  with  any  Tj.  Let  T  be  a  triangle  of  A.  Then 
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It  follows  that 


I  b  ♦ 

*-  rvtr  Tr 

p+q+r“4 


9(*+j) I 


■  Y  b  6 

^  W*Vir 


Hence  the  mapping  g  *-*■  b^  commutes  with  any  shift.  The  Lemma  is  proved. 


Corollary.  If  f  6  ij  and  X  e  A,  then  Xb^  is  invariant  under 
translates . 


Proof.  By  Lemma  5 


X<bIf(*+j)  -  bJf)  -  MTTj  -  If)  -  X(T(T^ f  -  f)) 


However,  T.f  -  f  e  w.j  hence  X(T(T.f  -  f))  ■  0.  This  shows  that 

j  4  j 


Now  let 


Xb 

If(‘+j)  - 

XbIf 

for  any 

j  e  *2  . 
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and 


N-1  N-1  4  3  . 
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m-0  n-0  i-1  j-1 
Since  £  ker  uh#  we  have 


ij"ij  I01/h 


with  h  >  0  and  N  -  |J/hJ  •  ( 


dist(f,Sh)  >  dist(f ,  ker  uh>  «  | w^f | / I I 


By  the  above  Corollary 
Jh 


“  1,2  l  l  AM'°)loi/hf  for  f  6  W« 

i-1  j-i  ,/n 


2  3  5 

For  f  i  x  H  xix2  w  have  °i/hf  *  h  f.  Hence 
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It  is  easy  to  verify  that 
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*  ‘bIf<!'i> 


bIf<-  W"  -  [bIf<M>  -  bH«-  }•!» 


Let  T  be  the  triangle  with  vertices  O  *  (0,0).  V  ■  (0,1)  and 


W  -  (1,1).  Then 

(x^,X2)  ■  u(0,0)  +  v(0,1)  +  w(1,1)  with  u+v+w  ®  1 
It  follows  that 

u  *  1-Xj  ,  v  ■  x2_x1  and  w  ■  x1 


Hence 
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and 
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Thus 


ijk  i  1  j  Ik!  V  12  1 
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r4-q,q-p,p 
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By  Lemma  5  we  have 

I,f,-*"V joJ.b«(5-5-,,Wb.P  • 

Therefore, 
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On  the  other  hand 
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Now  we  consider  another  triangle  with  vertices  U 

(1,0)  and  W  *  (1,1).  Then 


(0,0), 


u  •  1  -  x^,  v  ■  x^  “x2  an<^  w  “  x2  1 


x  «  (1,  Jli) 

Xi jk  l4'  4  ’ 
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Moreover,  we  have 


3  3  3  4 
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It  follows  that 
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In  conclusion  we  obtain 


4  3 
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a  X(0'0)If 
ij  ij 
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1 
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Thus 


lMhf  I  *  e  b^’2  ^  T2  h3  f°r  h  <  ’ 


Furthermore,  we  have,  for  any  g  e  C(IT), 
N-1  N-1  4  3 


u.,  -  I  l  l  i  ...  , 

”  Jio  nio  lit  ji,  11  ,/h 


N-1 

l  <*>, 


(m,-  7>-b  (m+  7,-  4)+b  <m+7»”  7)-bT_  (m+4*~ 
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1, 
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Nv  1 

+  l  C-bT„  (0,n)+b  (0,n+  4>-b_  (-  4,n+  7>+bTrt  (-  7,n+  f)] 
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(n-  4*n+  4)] 


It  is  easily  seen  that 


lbMJ  ‘  const  »  c 

l/h 


where  the  const  is  independent  of  h.  Hence  (17)  implies  that 


1 


lwhgl  *  4N  const  g  c  4  const  •  —  g  c 


This  shows  that 


“h  ‘  0(h’  • 


in 


Now  (14),  (16)  and  (18)  yield  the  desired  results 

4 


dist(f,S,  )  >  const  h 
h 


2  3 


for  some  positive  constant  and  the  function  f  »  xH  xix2 
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